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ABSTRACT 

We report on spectroscopic imaging observations of the nearby (z=0.022) galaxy cluster 
3C 129 performed with the ACIS detector on board of the Chandra X-ray observatory. Ap- 
plying a deprojection analysis which fully takes into account the spatially resolved X-ray energy 
spectra we investigate the radial variation of temperature and particle density of the Intracluster 
Medium (ICM). While earlier data indicated a weak cooling flow and a cool cluster core, the 
Chandra data suggests that the core is hotter than the rest of the cluster: at the center we infer 
7.9^2 4 keV while the mean temperature of all the other ICM shells is 5.1±0.3 keV (all errors 
on 90% confidence level). Based on the radial entropy profile we discuss the evidence for shock 
heating of the ICM at the cluster center. We discuss the possibility to use imaging spectroscopy 
data of more regular clusters to extract maps of the gravitational dark matter potential. 



Subject headings: galaxies: clusters: individual (3C 129) 
clusters 

1. Introduction 

While observations with the ROSAT and ASCA 
X-ray telescopes revealed non-homogeneous tem- 
perature of the ICM for a number of clusters (see 
e.g. Honda et al. 1996), the large scale radial 
temperature structure remained uncertain due to 
ROSAT's limited spectral resolution and cover- 
age, and ASCA's non-negligible energy dependent 
point spread function. Based on ASCA obser- 
vations of 30 galaxy clusters, Markevitch et al. 
(1998) found evidence for a general temperature 
decline with distance from cluster center. How- 
ever, independent ASCA studies of the same and 
other galaxy clusters did not confirm this result 
(e.g., White 2000, Kikuchi et al. 1999). The new 
X-ray observatories Chandra and XMM-Newton 
with broad energy coverage of 0.3-10 keV (Chan- 
dra) and 0.1-15 keV (XMM-Newton) and angular 
resolutions of 0.5 arcsec (Chandra) and 15 arc- 
sec (XMM-Newton) make it now possible to as- 
sess spatially resolved ICM X-ray energy spec- 
tra with sufficiently broad energy coverage and 
excellent angular resolutions. First observations 
of temperature or metallicity variations have al- 
ready been reported for a number of clusters, e.g. 
Coma (Arnaud et al. 2000), A 2142 (Markevitch 
et al. 2000), Hydra A (McNamara et al. 2000), 
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cosmology: dark matter — X-rays: galaxies: 

A 1795 (Tamura et al. 2001), Sersic 159-03 (Kaas- 
tra et al. 2001), and Abell 1835 (Schmidt et al. 
2001). On theoretical grounds one expects that 
the formation of clusters gives rise to substan- 
tial temperature, entropy, and metallicity gradi- 
ents. The interplay of gravitational collapse and 
heating by accretion shocks produces characteris- 
tic radial temperature and entropy profiles (e.g. 
Evrard 1990). Depending on when the ICM is en- 
riched with metals, also radial metallicity gradi- 
ents could be observed (Metzler & Evrard 1994). 
Furthermore, merger of major cluster components 
produce pronounced ICM non-uniformities, and 
hydro/N-body simulations predict that the gradi- 
ents could be observed several 10 9 years after the 
onset of the merger event (Roettiger et al. 1996). 
The decay of these non-uniformities might be re- 
duced by magnetic fields which organize the ICM 
into a filamentary structure and thereby substan- 
tially reduce the ICM's heat conductivity and par- 
ticle diffusion (Eilek 1999, but see also Narayan & 
Medvedev 2001). 

In this letter we describe Chandra observations 
of the nearby rich cluster of galaxies 3C 129. Due 
to its low galactic latitude (/=160°.5, 6=0.3°) the 
cluster has not been studied intensively at optical 
wavelengths. Based on ROSAT, Einstein, and EX- 
OSAT data, Leahy & Yin (2000) estimated a total 
0.2-10 keV luminosity of 2.7xl0 44 ergs s _1 , a total 
ICM gas mass of 3.6 xlO 13 M Q , and a total cluster 



mass of ^5 x 10 M©. The cluster harbors two 
radio galaxies: the prototypical head-tail galaxy 
3C 129, and the weaker FR I source 3C 129.1. 
The main emphasis of the Chandra observation 
had been a study of the pressure balance between S? 
the ICM and the radio plasma and this study will a 
be presented elsewhere. | 

In this letter we present a "spectroscopic °, 
deprojection-analysis" of the Chandra data. The 8 
deprojection-approach has been the standard- * 
technique to convert observed X-ray surface ^ 
brightness profiles into temperature and particle Q 
density profiles (see Fabian 1994, and references 
therein). Lacking spatially resolved X-ray energy 
spectra, these studies used several assumptions, 
namely that the ICM had a certain geometry, that 
the ICM was in hydrostatic equilibrium, and that 
the dark matter gravitational potential could be 
described by a certain shape and depth, together 
with an ICM temperature and metal abundance 
estimate derived from a non-imaging X-ray obser- 
vation. The spatially resolved energy spectra from 
the Chandra and XMM-Newton observatories now 
make it possible to directly derive the radial tem- 
perature, particle density and metallicity profiles 
based only on assumptions about the ICM geom- 
etry. A similar analysis has independently been 
developed by Allen et al. (2001). 

The rest of this letter is organized as follows. 
After describing the data set in Sect. 2 and the 
spectroscopic deprojection-method in Sect. 3, we 
will present the results of the observations in 
Sect. 4, and discuss their implications in Sect. 5. 
In the following we use Ho = 65 km s _1 Mpc -1 
and qo — 0.5; the cluster is thus a distance of 
100 Mpc and 1 arcmin corresponds to 28.7 kpc. 

2. Data Set and Data Preparation 

We use two 3C 129 pointings, one 30 ksec point- 
ing with the ACIS S CCDs taken on December 9th, 
2000, and one 10 ksec pointing with the ACIS I 
CCDs taken on January 9th, 2001. We detect 
X-ray emission near the radio cores of the radio 
galaxies 3C 129 and 3C 129.1 and exclude the cor- 
responding sky regions, as well as an additional 
point source in the south-east of the cluster cen- 
ter from the following analysis. The radio galaxy 
3C 129 has a 15 arcmin long tail which covers a 
substantial fraction of the field of view of the first 
observation at radial distances from 10 arcmin to 
25 arcmin from the cluster center. We exclude the 
corresponding sky region from the analysis. 

The analysis is performed with the CIAO 2 . 1 
software. We divide the 8.4' x 8.4' solid angle re- 
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Fig. 1— The ACIS I image (0.5-7.5 keV) with a 
32 x 32-pixel binning. The 4 ellipses are shown which 
have been used to characterize the symmetry of the 
ICM's surface brightness. The contour lines indicate 
the location of the radio galaxy 3C 129.1 at the cluster 
center and 3C 129 in the west of the cluster, and have 
been derived from 330 MHz VLA observations (Lane 
et al. 2001). The small rectangles cover the ACIS S 
and ACIS I field of views and show the sky-regions 
used for the deprojection analysis. The deprojection 
analysis uses only the regions with their centers within 
the outer ellipse. 

gions of the sky covered by each ACIS chip in 4x4 
"analysis regions" of each 126"xl26" (256x256 
pixels) solid angle coverage. This choice is a com- 
promise between keeping the analysis feasible and 
using detector response functions which are as lo- 
cal as possible. We construct background data 
sets with the tool make_acisbg using the blank 
field background data sets of M. Markevitch. Only 
events in the energy range from 1.2 keV to 7.5 keV 
are used. We use a rather high lower threshold 
to minimize the influence of uncertainties in the 
detector response and in the contribution of the 
diffuse galactic X-ray emission. Above 7.5 keV 
the signal to noise ratio is poor. Application of 
the same requirements on the background rate as 
used for the production of the background data 
sets, reduces the exposure time of the first and 
second pointing to 19.5 ksec and 9.3 ksec, respec- 
tively. We use the tools calcarf and calcrmf of 
A. Vikhlinin to compute for each analysis region 
an Auxiliary Response File (ARF) and a Response 
Matrix File (RMF) averaged over the correspond- 
ing detector area. The analysis uses the ACIS S 
chips S2-S4, and the ACIS I chips 10-14. Figure 
1 shows the binned ACIS I picture together with 
the ACIS S and ACIS I fields of view, the analysis 



regions, and the location of the two radio sources. 

3. Spectroscopic Deprojection Analysis 

The iso-brightness contours of 3C 129 show a 
slightly elliptical shape with the major axis along 
the east-west direction. Using the Chandra data 
together with the Einstein IPC contour map of 
(Leahy & Yin 2000), we define 4 ellipses with dif- 
ferent centers and elasticities to parameterize the 
ICM surface brightness distribution at different 
distances from the cluster center (see Fig. 1). 
Assuming rotational symmetry around the ma- 
jor axis of each ellipse and using an interpola- 
tion scheme we define 1 1 approximately ellipsoidal 
ICM shells. For each shell s (s =1, 2, 11) we 
compute the volumes V(s,r) along the line of sight 
of the 128 x 128 pixel analysis regions r (r =1, 
2, 95) by simple numerical integration. For 
each shell we define one Raymond-Smith plasma 
model described by one set of parameters tempera- 
ture, T(s), particle density, N(s), and metal abun- 
dances, A(s). The algorithm works from shell to 
shell inwards, fitting the plasma parameters of the 
shell under consideration and using for the outer 
shells the best fit parameters determined before. 
The expression minimized by the fit algorithm for 
shell s (counting from the cluster center outwards) 
reads: 

X '(T(.,),N(.,),A( s y,s)) = 

r'6s I' \ s'=s 

K (N K (r');E') j RS (T(s'), N(s>),A(s');E>) x 

D T ,(I',E')dE'^ (1) 

The first sum of the right term runs over the re- 
gions r' which have their centers within the shell 
s under consideration; the second sum runs over 
the energy bins I' of each energy spectrum; in the 
third sum s' runs over the shell under considera- 
tion s and all outer shells. C(r',I') is the number 
of counts found in energy bin /' of region r' . The 
constants c(r') ~ 1 allow for a slight deviation of 
the emissivity from the value expected from the as- 
sumed ellipsoidal geometry. We adjust these con- 
stants after fitting the flux normalization of a shell. 
The value £>l denotes the luminosity distance and 
k and j'rm describe the photo-electric absorption 
and the emissivity of the Raymond-Smith plasma 
as function of photon energy E', respectively. Fi- 
nally, D v i (/',£") is the detector response matrix 



which gives the detection area times the probabil- 
ity with which an X-ray photon of energy E' will 
be reconstructed within the energy bin I' (aver- 
aged over the detector coordinates of region r'). 
Due to the location of the cluster near the galac- 
tic plane we anticipated substantial variation of 
the neutral hydrogen column density across the 
field of view. We therefore performed radio obser- 
vations with the Dominion Astrophysical Obser- 
vatory to obtain ./Vjj-values with arcminute res- 
olution (principal investigator: T. Willis). Al- 
though the An-values change by 25% over the 
1.25° x2. 5° field of view, the Chandra observations 
lie in a region of rather constant values of between 
0.84xl0 22 cm" 2 and 0.9xl0 22 cnT 2 . We use in 
the following the -/Vu-values from the radio obser- 
vations. Independent determination of the mean 
iVR-value from Chandra data gives consistent re- 
sults. Note that our inferred hydrogen column 
densities deviate from those inferred by Leahy et 
al. (2000) {N H = (0.57±0.03)xl0 22 enr 2 ) from 
fitting the combined ROSAT and EXOSAT data 
with a one-temperature Raymond-Smith model. 
For each shell s we quote as radial distance R s the 
mean distances of all analysis regions contribut- 
ing to the plasma fit of this shell. We define the 
distance of a point from the cluster center as the 
half-diameter of the major axis of the associated 
ellipse. 

Technically the deprojection-analysis is per- 
formed with the CIAO-tool sherpa running on 
a single script which contains all the informa- 
tion about the geometry of the shells and anal- 
ysis regions. All errors are quoted on 90% con- 
fidence level. We estimate the increase of the 
statistical error on the parameters of the inner 
shells due to the statistical uncertainty of the fit- 
parameters of the outer shells with a simple Monte 
Carlo simulation and accordingly scale the errors 
on the fit-parameters. The uncertainty of the 
background normalization is estimated to be 10% 
(Vikhlinin ct al. 2001). Determining the associ- 
ated error in the parameters by repeating the full 
deprojection analysis with a background scaled up 
and down by 10%, we add these errors in quadra- 
ture to the statistical ones. The results of the 
outermost ICM shell depend on the assumed shell 
thickness. Since this is an inherent limitation of 
the deprojection method we do not show the fit re- 
sults of the 11th ICM shell in the following. The 
fits give satisfactory x 2 " va hies for 95 out of the 
104 analysis regions. For 9 regions large x 2_v alues 
indicate a deviation of the plasma properties from 
the ellipsoidal shell or contamination by a field 
source, and we exclude these regions from the over- 



all fit. 

4. Results 

The Chandra image shows statistically signifi- 
cant evidence for X-ray emission out to 19 arcmin 
from the cluster core. We do not find evidence for 
abrupt changes in the ICM surface brightness, so 
we have no direct evidence for shocks or contact 
discontinuities. 

Figure 2 shows the results of the deprojection 
analysis. Averaged over all shells, the mean clus- 
ter temperature is 5.1±0.3 keV. We have weak 
evidence for a hotter cluster core with k^T = 
7.9~^2A keV at R\ = 1.4 arcmin and a temper- 
ature decrease at the outermost 10th shell with 
fceT = 3.5jZQg keV Rio — 18.1 arcmin (fcs is 
the Boltzmann constant). Various estimates show 
that the low temperature of the outermost shell 
might be an artifact due to the soft diffuse galac- 
tic X-ray emission which begins to substantially 
contaminate the cluster emission at these dis- 
tances from the cluster center. The dotted line of 
Fig. 2(a) shows the projected temperatures deter- 
mined by fits of one component Raymond- Smith 
models. Naturally, the deprojected temperature 
varies more than the projected temperature which 
averages over several ICM shells. 

The particle density (Fig. 2(b)) decreases 
monotonically from 3.1 xlCP 3 cm~ 2 at the cluster 
center to 3.9 xlCP 4 cm~ 3 at the cluster periphery. 
The errors on the metal abundances of individual 
ICM shells are rather large (typically about 0.4 
solar abundances) and we do not show them here. 
The mean value averaged over all shells is 0.2±0.1 
solar abundance. The thermal gas pressure, de- 
fined as p = N k-gT is shown in Fig. 2(c). The 
entropy per particle relative to the mean value of 
all shells s is shown in Fig. 2(d) and has been 
computed from s = | fee In (Tp -2 / 3 ) - so with 
p « 0.6 mp TV and mp the proton mass. The 
entropy profile will be discussed in more detail 
in the next section. Finally, the radiative cool- 
ing times are shown in Fig. 2(e). The cooling 
times have been computed according to i coo i = 
| N k B T (A(iV,T))- 1 with the cooling func- 
tion A which denotes the plasma emissivity includ- 
ing thermal bremsstrahlung and line emission pro- 
cesses (Sutherland & Dopita 1993). All radiative 
cooling times exceed the Hubble time, the shortest 
cooling time being 30±g° Gyr. 

5. Discussion 

Based on a deprojection analysis of ROSAT sur- 
face brightness data Leahy & Yin (2000) inferred 
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Fig. 2. — Results of the deprojection analysis showing 
from top to bottom: (a) temperature T, (b) particle 
density N, (c) gas pressure, (d) entropy, and (e) ra- 
diative cooling times. The dotted line in panel (a) 
shows the projected temperature, and the dashed line 
in panel (d) is for illustrative purposes only (see text). 

a cluster temperature decreasing from ~6 keV at 
the cluster periphery (at ~20 arcmin from the 
cluster core) to about ~3 keV at the cluster cen- 
ter. The Chandra imaging spectroscopy data how- 
ever, does not confirm this temperature decrease. 
On the contrary, both, the projected and depro- 
jected temperature profiles indicate a hot clus- 
ter core. The discrepancy might be explained by 
one or more of the numerous assumptions entering 
the deprojection of surface brightness images men- 
tioned in the introduction, as e.g. a true dark mat- 
ter core radius substantially deviating from the as- 
sumed one. 

A hot cluster core might be the result of shock 
heating. This possibility can be studied with the 
entropy per particle which is able to distinguish 
between ICM heating by adiabatic compression or 
by shocks. By definition, adiabatic compression 
does not increase the ICM entropy, while shock 



heating, converting the energy of bulk plasma mo- 
tion into heat, does increase it. A possible inter- 
pretation of the entropy profile is based on the 
dashed line in Fig. 2(d) which shows a fit to the 
entropy of the second to tenth shell and suggests 
that the entropy increases with distance from the 
cluster center. Indeed, in the framework of hierar- 
chical clustering scenarios one expects that cluster 
formation through mass accretion produces a radi- 
ally increasing ICM entropy as consequence of the 
growing strength of the cluster accretion shock as 
the total virialized cluster mass accumulates (see 
e.g. Tozzi et al. 2001, and references therein). The 
large entropy of the cluster core might be the re- 
sult of shock heating by a jet or wind from the 
central radio galaxy 3C 129.1. The required power 
is modest: ~ 10 43 ^Zi ergs s" 1 where 

A(fce T) is the temperature increase at the clus- 
ter center and £30129.1 denotes the lifetime of the 
source. However, at the cluster center we do not 
find structure in the ICM surface brightness in di- 
rect support of this hypothesis. 

The analysis presented here is plagued by large 
statistical and systematic errors. The first stem 
from the modest integration time and the low 
surface brightness of 3C 129. The latter mainly 
results from the asymmetric shape of the ICM 
and possible contribution of the diffuse galactic 
X-ray emission at larger cluster core distances. 
For brighter, more symmetric clusters the spec- 
troscopic deprojection technique opens the possi- 
bility to determine the shape of the dark matter 
potential. The only uncertainties of such an anal- 
ysis derive from the uncertain ICM geometry and 
from undetected ICM pressure components as ki- 
netic pressure, magnetic field pressure and Cosmic 
Ray pressure. 
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